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Abstract. We present spectral and timing results from a long (130 ks) XMM-Newton EPIC observation of the 
nucleus of the Seyfert/LINER galaxy M81. During the observation the X-ray flux varied by 20%, but there was 
no significant change in spectral shape. The 2-10 keV spectrum is well described by a power law continuum and 
three narrow Fe Ka emission lines at 6.4, 6.7 and 6.96 keV. The three emission lines have equivalent widths of 
39^12, 47lj3 ^^'^ •^'^-15 respectively. The ratios of the three lines are thus more similar to those observed 
from the Galactic Centre region than to those typically observed from Seyfert galaxies. The high ionization lines 
most likely originate either from photoionized gas within 0.1 pc of the nucleus of M 81, or from a non-thermal 
distribution of cosmic-ray electrons interacting with the 0.2-0.6 keV thermal plasma which is found in the bulge 
of M 81. 
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1. Introduction 

At a distance of 3.6 Mpc (Freedman et al. 1994), the Sab 
spiral M81 is the closest galaxy showing the spectral sig- 
natures of a "low ionization nuclear emission-line region" 
(LINER). LINER characteristics are found in a signifi- 
cant fraction of all galaxies (between 1/5 and 1/3), and 
the nature of the relationship between LINERs and low- 
luminosity AGN (LLAGN) has been a long standing ques- 
tion (Ho, Filippenko & Sargent 1997, Heckman 1980). The 
presence of a LLAGN in M81 is well established: the un- 
resolved nuclear source emits a non-stellar UV continuum, 
broad optical/UV emission lines, and a significant flux 
of non-thermal X-rays (Ho, Filippenko & Sargent 1996). 
From dynamical studies, the mass of the central compact 
object is found to be w 4 x 10^ Mq (Ho 1999). With a 
bolometric luminosity of ~ 10**^ erg s~^, M81 is radiating 
at < 10^^ of its Eddington luminosity (Ho, Filippenko & 
Sargent 1996). 

The nuclear X-ray source has a power law spectral 
shape in the 0.3-10 keV energy range (Petre et al. 1993, 
Ishisaki et al. 1996, Swartz et al. 2003). The spectrum is 
absorbed by a small column density (~ 3 x 10^" cm~^) 
of material intrinsic to M81, and below 2 keV a multi- 
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tude of emission lines from circumnuclear, thermal gas 
can be seen in the XMM-Newton RGS spectrum of the 
nucleus (Page et al. 2003). Above 2 keV, the most remark- 
able feature seen in the spectrum is an emission line at 
6.7±0.1 keV, corresponding to Fe XXV Ka (Ishisaki et al. 
1996, Pellegrini et al. 2000). Using data from BeppoSAX, 
Pellegrini et al. (2000) claim that there is an absorption 
edge at ~ 8.6 keV, which could also arise in Fe XXV, 
indicating that the nucleus is seen through a large col- 
umn density (^ 2 x 10^^ cm~^) of highly ionized gas. 
The 2-10 keV luminosity is known to vary in the range 
(1.5 — 6) x 10'*" erg s"^ on a timescale of years, and vari- 
ations of about 30% have been observed on a timescale of 
10^ s (Pellegrini et al. 2000; Ishisaki et al. 1996). 

In this paper we examine the time variation and the 
2-10 keV spectrum of the M81 nucleus using the XMM- 
Newton European Photon Imaging Camera (EPIC) data 
from our long (130 ks) observation which was carried 
out as part of the RGS consortium guaranteed time pro- 
gramme. 

2. Observations and data analysis 

M81 was observed with XMM-Newton on 2001 April 22. 
The EPIC PN camera (Striider et al. 2001) was operated 
in small window mode for a single exposure of 130 ks du- 
ration. Each of the two EPIC MOS cameras (Turner et al. 
2001) were exposed for 122 ks, split into two exposures. 
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MOSl was operated in timing mode, while M0S2 was op- 
erated in full-frame imaging mode. AH three EPIC cam- 
eras were used in conjunction with medium filters. The 
nuclear X-ray source suffers from significant photon pile- 
up in the M0S2 data, and so these data are not included 
in our analysis. Pile-up is negligible for the nuclear source 
in the MOSI and PN exposures. 

The MOSl data were processed with SAS v4.1. For 
technical reasons, it was not possible to reprocess the 
PN data using SAS v4.1, so they were processed using 
SAS v5.3.3. In order to construct hghtcurves and spectra, 
source counts were extracted from a 44" radius circular 
region in the PN camera, and from a 39" wide strip in 
the MOSl camera. Both regions were centred on the nu- 
cleus and correspond to an encircled energy fraction of 
88% of the nuclear flux. In high resolution Chandra imag- 
ing of M81, the total X-ray flux originating from diffuse 
emission and point sources between 10" and 44" from the 
nucleus is only ~ 2% of the nuclear flux (Swartz et al. 
2003). Hence the EPIC spectra are completely dominated 
by emission arising from the nucleus and its immediate 
environment (within ~ 200pc). In both cameras, back- 
ground counts were taken from regions which were free 
from bright sources. In the PN camera these were on the 
same chip as the source region, and in the MOSl camera 
they were in the inner parts of the outer chips. 

Lightcurves were constructed from both PN and MOS 
1 data in the energy bands 0.3-2.0 keV and 2-10 keV 
with 1000 s time bins. To maximise the signal to noise, 
the lightcurves from the separate instruments were added 
to produce a single lightcurve for each energy band. 

When extracting spectra, we excluded periods of high 
background flux. This reduces the durations of the ex- 
posures to 119 ks for the PN camera and 93 ks for 
MOSl. Response matrices and effective area files specific 
to the source regions were constructed using the standard 
RMFGEN and ARFGEN SAS tasks. To optimise signal to noise 
the spectra from MOSl and PN were combined using the 
method of Page, Davis and Salvi (2003) in 45 eV bins 
which sample the EPIC spectral resolution well. The re- 
sultant spectrum has ^ 400 counts per bin at 6.4 keV 
and > 40 counts in every bin from 2 to 10 keV. Spectral 
analysis was carried out with XSPEC 11.1.0 (Arnaud 1996). 
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Fig. 1. EPIC lightcurves with 1000 s time bins in the 
0.3-2.0 keV band (top panel) and 2-10 keV band (mid- 
dle panel). The bottom panel shows the ratio of the two 
lightcurves; the dashed line shows a constant fit to the 
hardness ratio (see Section |OJ. 



X^/i/ = 141/125, corresponding to a null hypothesis prob- 
ability of 15%. This means that the low level variability is 
not accompanied by any significant gross spectral changes. 



3. Results 

3.1. Time variability 

The background-subtracted EPIC PN-l-MOSl hghtcurves 
for the 0.3-2.0 keV and 2-10 keV bands are shown in 
the top and middle panels of Fig. ^ respectively. The 
lightcurves show flickering, similar in character to the 
stochastic variability commonly seen in Seyfert galaxies 
(e.g. Green, McHardy & Lehto 1993), but of relatively 
low amphtude (~ 20% peak to trough). The ratio of the 
2-10 keV countrate to the 0.3-2 keV countrate is shown 
in the bottom panel of FigQ] A constant fit to these 
data points (shown as dashed line in Fignj results in 



3.2. Spectral study 

The 2-10 keV EPIC spectrum of M 81 is shown in Fig. El 
We began by fitting a power law model, absorbed by neu- 
tral material in our own Galaxy and within M 81, with the 
total column density fixed at the value of 7.6 x 10^° cm~^, 
as determined from the simultaneous RGS observation 
(Page et al. 2003). This model, and the data/model ra- 
tio are shown in Fig. |21 The parameters for this fit and 
all subsequent fits are given in Table While the overall 
continuum shape is well reproduced by the power law, the 
is very poor {x^/i' = 323/175), and the data show a 
significant excess over the model between 6.25 and 7.15 
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keV. The energy range of the excess emission corresponds 
to the Ka hnes of Fe I - Fe XXVI, although a photoelec- 
tric edge from Fe I at 7.1 keV could also contribute to the 
excess if there is significant reflection from cold material. 
Therefore we added a narrow 6.40 keV Gaussian emission 
line corresponding to Fe I Ka and a neutral Compton re- 
flection component (Magdziarz & Zdziarski 1995; PEXRAV 
in XSPEC) to the absorbed power law model. Both these 
components could arise in a cold reflecting medium, e.g. 
an accretion disc or a distant molecular torus (George & 
Fabian 1991; Matt, Perola & Piro 1991); the incHnation 
of the reflector was fixed at 60 deg to the observer's line 
of sight in the fit. The fit is much improved (A^^ of 69 
for 2 extra free parameters), compared to the absorbed 
power law model. However, this model is still rejected with 
> 99.9% confidence, and the data showed a significant ex- 
cess over the model at ~ 6.7 keV. Therefore we added a 
second narrow Gaussian emission line, at 6.68 keV, cor- 
responding to Fe XXV. The fit is improved significantly 
with the addition of the second emission line (Ax^ = 30 
for 1 extra free parameter). Although the model is now 
not rejected at the the 99% confidence level, it is still re- 
jected at the 95% level, and the data still show an excess 
over the model between 6.8 and 7.15 keV. Therefore we in- 
cluded a third narrow Gaussian emission line at 6.96 keV, 
corresponding to Fe XXVI, in the model. This results in 
another significant improvement in the fit (A^^ = 17 for 
1 extra free parameter) and the fit is now acceptable at 
the 95% confidence level. However, as shown in Tabled 
when the three Fe Ka emission lines are included in the 
fit, the amplitude of the reflection component is consistent 
with zero. We therefore removed the reflection component 
from the model and repeated the fit. The is the same 
for this model fit as it is when the reflection component 
is included, but with one fewer free parameter. Hence our 
best fit model is a power law with three Gaussian emission 
lines corresponding to Fe I, Fe XXV and Fe XXVI; this 
model is shown along with the data in Fig.O The 2-10 keV 
flux inferred from the fit is 1.1 x 10~^^ erg s~^ cm~^. To 
further check the consistency of the line associations with 
Fe I, Fe XXV and Fe XXVI we repeated the fit allowing 
the energies of the emission lines to vary. This results in 
a negligible improvement in goodness of fit (A^^ — 1 for 
3 extra parameters) and the emission line parameters are 
all consistent with those obtained when the line energies 
were fixed (see Fig.QJ. 

An alternative interpretation for the 6.25-7.15 keV 
emission is that it is produced by a single broad emis- 
sion line rather than by a combination of narrow lines. 
When a single broad Gaussian is tried, we obtain formally 
a good fit, with x^/i^ = 204/172. However, this model 
does not reproduce the sharp 6.25 keV shoulder of the 
emission feature, and is therefore less satisfactory than 
the multi-line fit (see Fig.O. The best fit central energy 
of the line is 6.69lo;o7 keV consistent with Fe XXV Ka, 
and the full width half maximum (FWHM) of the line is 
28000 ± 7000 kms-i. 




4 6 8 10 

Channel energy (keV) 



Fig. 2. Top panel: EPIC spectrum of M81 (datapoints) 
and the best fit power law model (stepped line). Bottom 
panel: ratio of the datapoints to the power law model. 



Pellegrini et al. (2000) claimed to have detected a pho- 
toelectric edge at 8.6lQ g keV with optical depth r = 
0.15^Q Qg in the spectrum of M81, which they hypothe- 
sized could be due to Fe XXV. We have searched for Fe K 
photoelectric edges in our spectrum by adding an edge to 
our best fit model. The threshold energy of the edge was 
allowed to vary between 7 and 9.3 keV, corresponding to 
K edges from Fe I - Fe XXVI. However, we find that the 

is unchanged by this extra model component, and the 
best fit optical depth for the edge is zero, regardless of 
its energy. To determine upper limits we stepped the edge 
energy and optical depth over a grid of values, to produce 
the confidence contour shown in Fig. [S| 

4. Discussion 

In common with previous investigations, we find that the 
2-10 keV spectrum is dominated by a power law contin- 
uum. The spectral slope we obtain (F = 1.88 ± 0.02) is 
consistent with that seen below 2 keV in the simultane- 
ous RGS observation (F = 1.94 ± 0.06; Page et al. 2003). 
Low level variability, of around 20% peak to trough is ob- 
served during the observation, but there are no detectable 
spectral changes. Similar fractional variations in inten- 
sity, over 1-2 day timescales, have also been observed in 
ASCA (Ishisaki et al. 1996) and BeppoSAX (Pellegrini et 
al. 2000) observations. The mean EPIC countrate during 
the XMM-Newton observation corresponds to a 2-10 keV 
flux of 1.1 X 10^^^ erg s~^ cm~^, around 30% of the flux 
at the time of the BeppoSAX observation. Thus the X- 
ray variability of M81 observed on a 1-2 day timescale is 
consistent with the linear relationship between flux and 
variability amplitude that is observed in Seyfert galax- 
ies, whereby the variability amplitude is proportional to 
the mean flux, and so the fractional variations change 
little from observation to observation even though the 
mean flux may change considerably (Uttley & McHardy 
2001). On the other hand, the variability amplitude of 
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Table 1. Results of spectral fits to the EPIC data. All fits were performed over the 2-10 keV energy range. All the 
spectral models included absorption from neutral material with a column density of 7.6 x 10^" cm^^. Uncertainties 
are quoted at 95% for 1 interesting parameter (i.e. Ax^ = 4). Apj^ is the power law normalisation in units of 
10~^ photons cm~^ s~^ keV~^. R is the solid angle of the sky covered by the reflector as seen by the X-ray source 
in units of 27r steradians. When a parameter was constrained by the limits of the fit-range rather than by A^^ the 
corresponding limit has been marked with an 



6.40 keV 6.68 keV 6.96 keV 

Model r Ap.L. E.W. E.W. E.W. R ^ jv Prob 

(eV) (eV) (eV) 

to 07 _ _ - _ 323/175 7.4 x 10"" 

-o'w 26ti^ - - l-Stae 245/173 2.5 x 10"'' 

iaw 30ti^ 38ti^ - l-Slo'e 218/172 1.0 x 10"^ 

tais 35^13 43ti^ 32^^^ 0.6j:[]:^^ 202/171 5.3 x 10"^ 

to'oT 39li2 47^13 37^15 - 202/172 5.8 x 10"^ 

- 160l4^ - - 204/172 4.8 x 10"^ 



P.L. 


J^-»D-0.02 


3.55: 


P.L.-HGau-l-refl 


1 05+0.05 


3.87: 


P.L.-|-2xGau+refl 


1 02+005 
^■^^-0.05 


3.78: 


P.L.-|-3xGau+rofl 


1 S8+0-05 


3.68: 


P.L.-|-3xGau 


1 ss+0.02 


3.64: 


P.L.+broad Gau^ 


1 SS+0-02 


3.65: 



f In this model the central energy and width a of the emission line are free parameters. The best fit values are 6.691g 

0.06 



and 0.27tn 06 keV respectively. 



M81 (and LINERs in general) is much lower than would 
be expected from extrapolating Seyfert variability ampli- 
tudes to lower luminosity (Ptak et al. 1998). The spectral 
index obtained from the XMM-Newton observation is con- 
sistent with those obtained from observations with ASCA 
(r = 1.85 ±0.04, Ishisaki et al. 1996) and with BeppoSAX 
(r = 1.86 ± 0.03, Pellegrini et al. 2000). Thus the contin- 
uum shape appears to be unaffected by intensity variations 
either over timescales of days or over timcscales of years. 

4.1. Emission and absorption from ionized Fe 

The high quality of our EPIC spectrum allows us to re- 
solve the Fe Ka line emission, which can be modelled ei- 
ther by a single broad Fe XXV emission line or as the 
sum of three narrow components: Fe XXVI, Fe XXV, and 
a component representing emission from lower ionisation 
Fe species. Fe Ka emission is commonly seen in Seyfert 
galaxies, and is thought to be due to reprocessing of the 
primary X-rays by surrounding material. If the material 
lies along the line of sight, we will see both absorption and 
emission features in the X-ray spectrum. This is the hy- 
pothesis favoured by Pellegrini et al. (2000), because their 
BeppoSAX spectrum appeared to show both an emission 
line and an absorption edge from Fe XXV. However, we 
can rule out a Fe XXV K edge with significant optical- 
depth (t > 0.1) at the > 99.7% confidence level in our 
EPIC spectrum. One possible explanation for this discrep- 
ancy might be that the X-ray flux at the time of the EPIC 
observation was only about 30% of the flux at the time of 
the BeppoSAX observation: assuming the absorber is pho- 
toionized, the corresponding reduction in the ionization 
level of the absorber might explain the lack of an Fe XXV 
edge in the EPIC spectrum. If this were the case, we would 
expect the dominant ionic species of Fe in the absorber to 



change from Fe XXV during the BeppoSAX observation 
to Fe XX - Fe XXIV in the XMM-Newton observation 
(Kallman & McCray 1982). Instead of the r ~ 0.2 Fe XXV 
K edge at 8.8 keV, there should be 1-4 Fe K edges between 
8.0 and 8.6 keV. The relative ionic abundances of Fe at this 
ionization level (Kallman & McCray 1982), and the cor- 
responding K-shell photoionization cross sections (Verner 
& Yakovlev 1995), imply that at least one of these should 
have T > 0.06 at threshold, but this is ruled out with 95% 
confidence by our upper limits given in Fig. [S] 

An absorber in which the dominant species of Fe are 
Fe XX - Fe XXIV should also produce a host of Fe L 
absorption lines in the soft X-ray. The deepest of these 
lines lie at 12.82 A, 12.29 A, 11.71 A, 11.00 A, and 10.64A 
for Fe XX - Fe XXIV respectively (e.g. Verner, Verner 
& Ferland 1996). If the column density implied by the 
Fe XXV K edge found by Pellegrini et al. (2000) were 
shared evenly between these ions, and assuming a veloc- 
ity dispersion of cr > 50 km s^^ in the absorber, we would 
expect these lines to have equivalent widths of at least 
29 mA, 27 mA, 19 mA, 15 mA and 13 mA respectively. We 
have looked for these absorption lines in the simultane- 
ous RGS spectrum of M 81 studied by Page et al. (2003) 
by repeating their spectral fit, adding these minimum ab- 
sorption lines one at a time to their best fit model, and 
measuring the Ax^- For every ion except Fe XXII the 
absorption line makes the poorer, by Ax^ values in 
the range 1.4 to 9. For Fe XXII, adding an absorption 
line of 19 mA equivalent width improves the x"^ by 0.9. 
Increasing the equivalent widths of any of these lines (in- 
cluding Fe XXII) from their expected minima results in 
an increased x^- Furthermore, in this range of ionization, 
3-20% of O should be in the form of O VIII. Assuming the 
absorber has solar abundances (Anders & Grevesse 1989), 
and a velocity distribution with tr > 50 km s~^, the ab- 
sorber should give rise to O VIII Lya and Ly/3 absorption 
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Fig. 3. Top panel: EPIC spectrum of M81 (datapoints), 
the best fit power law + three Fe emission lines model 
(black stepped line), and the power law + broad Gaussian 
model (grey stepped line). Middle panel: contribution of 
the datapoints to the multiplied by the sign (model- 
data) for the power law + three Fe emission lines model. 
Bottom panel: contribution of the datapoints to the 
multiplied by the sign (model-data) for the power law + 
broad Gaussian model. 



lines with equivalent widths of > 28 mA and > 16 mA re- 
spectively. Such absorption lines are not seen in the simul- 
taneous RGS spectrum (Page et al. 2003). Instead, O VIII 
Lya is seen in emission, and while this might mask the 
absorption line, the corresponding Ly/3 emission is much 
weaker. To see whether such a Ly/3 absorption line might 
be present in the RGS spectrum, we have again repeated 
the spectral fit of Page et al. (2003) to the RGS spectrum, 
adding an absorption line at the position of O VIII Ly/3 to 
their best fit model. The best fit equivalent width is zero, 
and an absorption line of 16 mA equivalent width results 
in a Ax^ of 4, and so is excluded at 95% confidence. In fact 
we do not see any evidence for ionized absorption in the 
RGS spectrum. Therefore we conclude that the lack of an 
Fe XXV K absorption edge in our EPIC spectrum is not 
explained by the reduction of the ionization. Either the 
ionized absorber has moved out of the line of sight, or the 
2.8cr edge feature in the BeppoSAX spectrum was a sta- 
tistical fluctuation rather than a real feature. Therefore, if 
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Fig. 4. Confidence intervals for the energies and normal- 
isations of the three narrow emission lines used to fit the 
bump between 6.25 and 7.15 keV. The solid, dashed and 
dotted contours correspond to 1, 2 and 3(T confidence for 
two interesting parameters respectively. The solid vertical 
lines indicate the expected energy for Fe I, Fe XXV and 
Fe XXVI. 
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Fig. 5. Confidence contour (essentially upper limits) for 
a photoelectric edge due to ionized Fe. The solid, dashed, 
and dotted lines correspond to confidence intervals of 1, 2 
and 3cr respectively for two interesting parameters. 



the Fe Ka line(s) are reprocessed emission from the X-ray 
radiation, they come from outside the line of sight. 

If a broad Fe XXV line is responsible for all the 
emission, the large velocity width (FWHM = 28000 ± 
7000 km s~^) implies that the emission originates within 
a few hundred Schwarzschild radii of the compact object. 
Furthermore, because the emission feature has a relatively 
high equivalent width (160 ± 40 eV), it must have a sig- 
nificant Thomson optical depth (i.e. > 0.1) and cover a 
significant fraction of the sky as seen by the X-ray source 
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(see Bianchi & Matt 2002) . The inner part of an accretion 
disc is the only viable location for gas with this combi- 
nation of properties, and reflection from an accretion disc 
surface is widely hypothesized for broad Fe Ka emission in 
AGN. However, the inner part of an accretion disc is un- 
likely to be responsible for broad 6.7 keV line emission in 
M 81, on both theoretical and empirical grounds. Firstly, 
because of its low bolometric luminosity (L < 10^'^iEdd), 
M 81 probably contains an advective flow in the central re- 
gions rather than an accretion disc (Narayan & Yi 1995). 
This hypothesis is further supported by the lack of any 
soft, blackbody-like component in the soft X-ray spec- 
trum of M 81 (Page et al. 2003); such a component is 
commonly observed in Seyfert galaxies and is thought to 
be the high energy tail of the emission from the inner ac- 
cretion disc (Arnaud et al. 1985). The second reason that 
a broad Fe XXV line is unlikely to originate in an accre- 
tion disc in M 81 is that the X-ray spectrum lacks any of 
the other features expected from such a disc. In particular, 
the soft X-ray spectrum does not contain the strong, broad 
emission features predicted to occur in discs which are suf- 
ficiently ionized to produce significant Fe XXV emission 
(e.g. O VIII and O VII emission lines and radiative re- 
combination continua; see Nayakshin & Kallman 2001 and 
Ballantyne, Ross & Fabian 2001). Instead, the only emis- 
sion features observed in the soft X-ray spectrum arc the 
narrow emission lines from a spatially extended thermal 
component (Page et al. 2003). Thus a single broad Fe XXV 
emission line, as well as failing to reproduce the sharp 
6.25 keV shoulder of the emission feature, is physically 
inconsistent with the other X-ray and multiwavelength 
properties of M 81. However, for a novel explanation for 
a broad 6.7 keV emission line, arising at the interface 
between the thin, outer disc and the inner, radiatively- 
inefficient, hot flow, see Dewangan et al. (2004). 

Because the broad line model for the Fe Ka emis- 
sion appears unlikely on physical grounds, we consider the 
combination of narrow 6.4, 6.7 and 6.9 keV emission lines. 
The three components each contribute about a third of 
the line emission, and the combined equivalent width is 
120 ±20 eV. For the 6.4 keV emission line, which arises in 
cool material, the observed equivalent width of ~ 40 eV 
is rather low compared to the line width observed in most 
Seyfert galaxies (typically ~ 100 eV, e.g. Nandra et al. 
1997). If the line arises in a Compton thick, plane paral- 
lel reflector, then the small equivalent width implies either 
that the reflector is highly inclined to our line of sight (i.e. 
i > 80°) or that the reflector covers less than 50% of the 
sky as seen by the X-ray source (George & Fabian 1991). If 
the reflector has a funnel shape (e.g. the molecular torus of 
AGN uniflcation schemes) then the low equivalent width 
implies an even lower sky-coverage of the reflector (Matt 
et al. 1991). Alternatively the low equivalent width of the 
6.4 keV line could be explained if the reprocessing region 
is Compton-thin (Leahy & Crcighton 1993). In any case, 
the low equivalent width of the 6.4 keV line implies that 
in M81 there is not as much cold. X-ray illuminated ma- 
terial as there is in a typical Seyfert galaxy. In contrast, 



the 6.7 and 6.96 keV lines have a larger equivalent width 
in M81 than is normally observed in Seyfert 1 galaxies, 
which typically have 6.4 keV lines that are at least four 
times stronger than either the 6.7 or 6.96 keV line; in- 
deed the 6.7 or 6.96 keV lines are often not detectable in 
EPIC X-ray spectra of Seyfert Is (e.g. Pounds et al. 2003, 
Blustin et al. 2003). The calculations by Bianchi & Matt 
(2002) allow us to place some interesting constraints on 
photoionized material that might give rise to these lines. 
According to these calculations the column density of illu- 
minated, reflecting material must be > 10^^ cm^^ in order 
to produce the observed equivalent widths. To maintain a 
sufficiently high level of ionization in the illuminated ma- 
terial to produce the 6.7 and 6.96 keV lines requires that 
the material lies quite close to the illuminating source. 
From equation 6 of Bianchi & Matt (2002), and for the 
observed 2-10 keV luminosity of M 81 (1.7 x 10''° erg s^^), 
we find that the illuminated material responsible for the 
6.7 keV line must lie within 0.1 pc of the nucleus, while the 
material responsible for the 6.96 keV line must lie within 
0.01 pc. 

Although the 6.4-7 keV complex could not be sepa- 
rated into its components in the ASCA and BeppoSAX ob- 
servations, the total equivalent width measured with these 
data was consistent with our EPIC value, at 170 ± 60 eV 
(Ishisaki et al. 1996) and 95j:4^ eV (Pellegrini et al. 2000) 
respectively. Unfortunately, the uncertainties on these 
measurements are large enough that they are also con- 
sistent with constant flux from the lines (i.e. a reduction 
in the equivalent width by a factor of three between the 
BeppoSAX and XMM-Newton observations). 

4.2. Comparison to the Galactic Centre 

It is interesting to note that the Fe Ka line ratios are 
relatively similar to those observed in the 1 deg x 1 deg 
(150 pc X 150 pc) region around our own Galactic Centre 
(Tanaka et al. 2000). In particular, emission lines at 6.4 
keV, 6.7 keV and 6.9 keV are seen in both the Galactic 
Centre and in M81, and in both cases the 6.7 keV line is 
the strongest (Tableland Tanaka 2002). In the case of 
the Galactic Centre, the origin of the Fe K lines is contro- 
versial (e.g. Ebisawa et al. 2003; Predehl et al. 2003). The 
6.4 keV emission is thought to arise through fluorescence 
in cold Fe, which is either illuminated by X-rays emitted 
during a past, comparatively more active period of Sgr A* 
(e.g. Murakami, Koyama & Maeda 2001), or bombarded 
by low energy cosmic ray electrons (Valinia et al. 2000). 

Possible models for the 6.7 keV and 6.96 keV lines in- 
clude a high temperature {kT ~10 keV) thermal plasma 
(Koyama et al. 1989), charge exchange between cosmic- 
ray Fe ions and interstellar H (Tanaka, Miyaji & Hasinger 
1999), interaction of non-thermal electrons with a low 
temperature 0.3 keV) thermal plasma (Masai et al. 
2002; Dogiel et al. 2002; Valinia et al. 2000), or a large 
population of faint accreting binaries (Wang et al. 2002). 
However, there now appear to be signiflcant problems with 
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the first two of these models for the 6.7 and 6.96 keV 
Unes: a 10 keV thermal plasma is too hot to be gravi- 
tationally bound to the Galaxy, and so strong magnetic 
fields (~ 30/iG) would be required to confine the plasma 
(Tanuma et al. 1999), while charge exchange between cos- 
mic ray Fe ions and interstellar H is expected to produce 
8.5 and 8.9 keV lines rather than the 6.7 and 6.96 lines 
which are observed (Masai et al. 2002). 

Although the integrated contribution of accreting bi- 
nary stars may plausibly explain the 6.7 and 6.96 keV 
lines from the Galactic Centre region (Wang et al 2002), 
this would seem unlikely for M81. The population of dis- 
crete X-ray sources near the Galactic Centre has a ra- 
dial surface density which declines as where R is dis- 
tance from the centre (Muno ct al. 2003). This is sim- 
ilar to the spatial distribution of stars selected in the 
infrared (Serabyn & Morris 1996), suggesting that the 
X-ray source population is intimately linked to the un- 
derlying population of bulge stars. Furthermore, the lu- 
minosity function of the discrete sources is steep, with 
dN/dL oc L~^-'' , at the faintest fluxes probed by Muno 
et al. (2003) so that the emission from discrete sources 
is dominated by low luminosity objects (2-8 keV lumi- 
nosity < 10^^ erg s~^), likely to be low mass accreting 
binaries such as cataclysmic variables. Such objects show 
strong Fe XXV and Fe XXVI K a Unes (see e.g. Mukai 
et al. 2003; Wu, Cropper & Ramsay 2001) and will pre- 
dominantly be drawn from the old stellar bulge population 
because they are long lived (Warner 1995). The bulges of 
the Milky Way and M 81 have similar overall luminosities: 
5 X 10^ L0 for M81 (determined from V-band photome- 
try: Tcnjcs, Haud & Einasto 1998), and 6 - 11 x 10^ L© 
for the Milky Way (determined from near-IR photome- 
try: Kent, Dame & Fazio 1991; Freudenreich 1998). We 
would therefore expect similar integrated Fe K line lumi- 
nosities from the low luminosity binary populations in the 
central regions of M81 and the Milky Way. The luminos- 
ity of the Fe Ka lines emitted near the Galactic Centre 
is 1.1 ± 0.1 X 10^^ erg s^^ (Yamauclii ct al. 1990), and 
assuming a distance of 3.6 Mpc to M81, the Fc Ka line 
emission from M81 is 1.6 ± 0.3 x 10'^* erg s^^, about 150 
times more luminous than that from the Galactic Centre. 
Therefore, whatever fraction of the Fe K line emission in 
the Galactic Centre comes from accreting binaries, such 
sources are unlikely to make a substantial contribution to 
the Fe K lines in the central part of M81. 

On the other hand, if the Galactic Centre Fe lines 
originate from a non-thermal distribution of cosmic-ray 
electrons, interacting with both cold and hot (0.2-0.6 
keV) gas, the Fe lines in M81 could quite plausibly arise 
through the same mechanism. In the radio band it has 
already been noted that the Galactic Centre and M81 
appear to be remarkably similar sources, with M81 re- 
sembling a scaled-up version of the Galactic Centre, from 
the core-jet structures (M81* and Sgr A*) seen on sub-pc 
scales (Bietenholz, Bartel & Rupen 2000; Lo et al. 1998) 
to the large arc-like structures seen on tens to hundreds 
of pc scales (Kaufman et al. 1996; Yusef-Zadeh, Morris & 



Chance 1984). The radio luminosity of M81* is about lO'' 
times higher than that of Sgr A*, which suggests that the 
non-thermal electron population in the central regions of 
M 81 could be supplied with sufficient energy to power the 
Fe Ka emission. Furthermore, myriad soft X-ray emission 
lines from 0.2-0.6 keV thermal plasma have already been 
detected from the central kpc of M 81 where the most im- 
portant radio structures are found. 

The spatial extent of the 6.4, 6.7 and 6.9 kcV line emis- 
sion in M 81 will be a key diagnostic to determine whether 
they are analogous to those observed from the Galactic 
Centre or arc due to photoionization by the nuclear X- 
ray source. The lines in the Galactic Centre come from 
a region more than 100 pc across; if the M 81 lines are 
a scaled-up version of the same phenomenon, then they 
can be expected to arise from region of similar size or 
larger (i.e. extending out to at least 6 arcseconds from 
the nucleus of M81). On the other hand, if the lines are 
photoionized by the nucleus they must originate within 
0.1 pc of the nucleus (6 milli- arcseconds at the distance 
of M81). Unfortunately, at the spectral resolution of the 
EPIC CCDs the peak intensity of the lines is at most 
~ 20% of the underlying AGN continuum, and XMM- 
Newton^s spatial resolution is not sufficiently good that we 
can distinguish such a weak component of extended emis- 
sion. Chandra has much better spatial resolution, but it 
has a much lower sensitivity than EPIC around Fe K, and 
the existing observations suffer badly from photon pile-up 
close to the nuclear source. Currently, Fe K emission has 
been detected neither from the nuclear source, nor from 
the circumnuclear 10-30 arcsecond region (Swartz et al. 
2003). In principle, a very long observation using the high- 
energy transmission gratings on Chandra might allow the 
measurement of the spatial extent of the Fe K line emit- 
ting region in M81 because it could achieve the required 
spatial and spectral resolution without significant photon 
pile-up. 

5. Conclusions 

The 2-10 keV X-ray spectrum of the nucleus of M81 con- 
sists of a power law continuum and three Fe Ka emission 
lines at 6.4, 6.7 and 6.96 keV. The power law component 
has a photon index F ^ 1.9, consistent with the contin- 
uum observed below 2 keV with the RGS. The X-ray flux 
varied by 20% during the 130 ks observation but no gross 
spectral variations were detected. The spectrum does not 
require a reflection component, and shows no evidence for 
absorption edges due to highly ionized Fe. The equiva- 
lent widths of the 6.4, 6.7 and 6.96 keV emission lines 
are SdtH, 47^13 and 37^11 eV respectively. These equiv- 
alent widths are signiflcantly different to those normally 
observed in Seyfert galaxies, which typically have stronger 
6.4 keV lines and weaker 6.7 and 6.96 keV lines. However, 
the ratios of the three lines in M 81 are relatively simi- 
lar to the those observed from the Galactic Centre region. 
The nucleus of M 81 does not appear to contain as much 
cold, X-ray illuminated material as most Seyfert galax- 
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ies. The 6.7 and 6.96 keV emission lines might originate 
in photoionized gas within 0.1 pc of the nuclear contin- 
uum source. Alternatively, these emission lines might have 
a similar origin to the Galactic Centre Fe line emission, 
possibly through the interaction of non-thermal cosmic ray 
electrons with the 0.2-0.6 keV thermal plasma component 
of the M81 bulge. 
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